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ABSTRACT 


In  this  report,  we  discuss  the  amplitude,  location,  and 
the  relative  separation  estimation  accuracy  of  radar  targets 
closely  spaced  in  range  using  the  Cramer-Rao  bound.  It  is  assumed 
that  the  phases  of  successive  signals  are  coherent  and  therefore 
contain  relative  line-of -sight  location  information.  It  is  shown 
that  this  information  can  substantially  reduce  the  estimation 
error  when  compared  with  the  case  where  the  relative  signal 
phase  is  random. 
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1.  INTRODUCTION 


Resolution  of  closely  spaced  signals  has  been  a  subject  of 
interest  for  many  years,  see  for  example  [1]  -  [14].*  The 
traditional  signal  resolution  problem  deals  with  radar  signals, 

[1]  -  [10],  while  more  recent  emphasis  has  been  on  optical  signals 
[10]  -  [14],  The  main  difference  between  the  radar  and  optical 
signal  is  that  the  radar  signal  is  a  complex-valued  process  cor¬ 
rupted  by  complex-valued  Gaussian  white  noise  while  the  optical 
signal  is  a  real-valued  process  corrupted  by  a  real-valued  non- 
Gaussian  process  and  it  is  usually  a  random  signal  in  colored 
noise  problem. 

There  are  two  basic  issues  pertaining  to  the  problem  of 
signal  resolution,  i.e.,  the  recognition  of  the  existence  of 
multiple  signals  (the  problem  of  detection)  and  the  estimation 
o l  pertinent  signal  parameters  (the  problem  of  estimation)  .  The 
problem  which  has  drawn  more  attention  in  the  literature  is  the 
problem  of  detection.  This  is  because  the  ordinary  multiple 
:nr.!';'  likelihood  ratio  test  procedure  is  ill-conditioned  for  this 
}■■!• 'i»lc.*ai  (see  [1]  and  [3]),  and  alternate  information  criteria 
must  bo  employed  (see  [1]  -  [5]  and  [14]). 


m  list  is  not  intended  to  be  exhaustive.  Ref.  [1]  has  a 
i'  rt  eut  concise  review  of  the  open  literature. 


As  to  the  problem  of  estimation,  the  maximum  likelihood  es¬ 
timator  has  been  widely  used  due  to  its  asymptotic  optimality 
property  [15],  [16]  and  intimate  tie  with  the  detection 

algorithm.  [1]  -  [10].  The  difficulty  of  applying  the  maximum 
likelihood  estimator  is  in  the  computational  burden  since  it 
is  generally  a  multi-dimensional  optimization  problem. 

Furthermore,  there  may  exist  many  local  maxima  in  the  likelihood 
function  for  a  certain  class  of  signals.  Once  a  sufficiently 
close  initial  guess  is  found  however,  an  iterative  optimization 
algorithm  (e.g.,  the  gradient  algorithm)  can  usually  find  the 
optimal  solution. 

Another  difficulty  in  the  problem  of  estimation  is 
performance  evaluation.  Due  to  the  nonlinear  nature  of  the  pro¬ 
blem  an  exact  expression  for  the  error  covariance  of  the  estimates 
cannot  be  obtained.  Fortunately,  the  error  covariance  of  the 
maximum  likelihood  estimates  can  be  closely  predicted  using 
the  Cramer-Rao  bound  when  the  signal-to-noi se  ratio  is  large. 
Furthermore,  the  maximum  likelihood  estimates  are  asymptotically 
Gaussian  and  approach  the  Cramer-Rao  bound,  [15]  -  [16]  . 

For  these  reasons,  Cramer-Rao  bound  analysis  has  gained 
significant  popularity  in  recent  years  for  predicting  the  es¬ 
timation  performance  in  the  signal  resolution  problem,  [7]  -  [10], 
[12]  -  [13]  . 
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The  Cramer-Rao  bound  for  parameter  estimation  of  closely 
spaced  radar  signals  was  presented  in  [8]  -  [10].  Refs.  [8]  and 
[9]  used  some  simplifying  assumptions  in  which  the  amplitudes 
(and  phases  for  [8])  of  the  signals  were  assumed  known. 

Closed  form  expressions  were  obtained  for  the  two-signal  case 
which  provided  intuitive  insights  about  the  estimation  accuracy. 
These  bounds  are  unfortunately  too  optimistic.  In  Ref.  [10],  the 
Cramer-Rao  bound  for  jointly  estimating  amplitudes,  phases,  and 
time  delays  was  obtained.  This  provided  a  realistic  bound 
for  this  problem  and  the  close  agreement  with  simulation  results 
was  also  demonstrated.  In  this  report,  we  extend  the  work  of  [10] 
to  include  the  case  where  the  phase  angles  of  successive  signal 
are  coherent.  There  are  many  applications  for  which  this 
condition  is  satisfied.  For  example,  if  one  wishes  to  estimate 
the  locations  of  scattering  centers  on  a  distributed  target,  the 
phase  difference  of  two  successive  scatterers  is  proportional  to 
their  separation  along  the  line  of  sight.  Even  if  two  point 
targets  are  completely  resolved  in  terms  of  base  band  pulse  shape, 
their  relative  separation  is  still  contained  in  their  phase 
difference  given  that  the  radar  can  process  the  returned  signal 
coherently.  This  second  example  is  indeed  the  basis  of  the 
phase-derived  range  technique  illustrated  in  [18].  As  will  be 
’.awn  later,  this  phase  information  is  able  to  substantially 


reduce  the  estimation  error. 


T 


This  report  is  organized  as  follows.  In  the  next  section, 
we  introduce  the  signal  model  and  the  Cramer-Rao  bound.  Equations 
for  computing  the  Fisher  information  matrix  and  some  numerical 
results  are  presented  in  Section  3.  A  summary  and  conclusions 
are  stated  in  the  last  section.  Three  appendices  are  attached. 

In  Appendix  A,  we  discuss  the  Cramer-Rao  bound  for  multiple  radar 
pulses.  In  Appendix  B,  we  give  the  formulation  and  structure  of 
the  maximum  likelihood  estimator  for  resolving  closely  spaced 
targets.  In  Appendix  C,  we  re-state  the  Fisher  information 
matrix  of  Ref.  [10]  for  the  random  phase  angle  case  for  the  pur¬ 
pose  of  comparison. 
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2.  SIGNAL  MODEL  AND  THE  CRAMER -RAO  BOUND 

Let  s(t)  denote  the  complex  low-pass  transmitted  waveform. 
The  reflection  from  the  ith  point  target  (or  scatterer)  is 
a^  s(t-r^).  The  total  received  signal  from  a  group  of  N  targets 
may  be  expressed  by 

N 

r(t)  =  2  a.  s  ( t— t i )  +  n (t)  (2.1) 

i=l 

T  A 

where  (t^,...,t  )  =  t_  =  locations  of  point  targets  in  terms 

1  of  relative  range 

T  A 

(a^,...,a^)  =  a  =  complex  scattering  amplitudes 

n(t)  =  receiver  white  noise  which  has  a 

complex  Gaussian  distribution  with 
variance  Nq/2. 

The  log  likelihood  ratio  of  signals  present  versus  signals 
absent  is 

-  N 

InA  =  jji  / 1  r  (t)  -  £  a.s (t-i . ) ] 2dt  +  /|r(t)|2dt  (2.2) 

o  J  i=l  o  J 


dt 


dt 


bet  y  ^  =  J s ( t— t  ^ ) r* (t) 

pi j  =  /  S  (t-Ti)  S  (t-Tj ) 

A  r  ,T 

P  4  lpij] 

where  the  superscript  denotes  the  complex  conjugate  operation. 
The  log  likelihood  ratio  may  be  rewritten  as 


IT  T  T 

InA  =  (a/y  +  a*  y*  -  aTa*) 

o 


(2.3) 
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The  above  constitutes  the  signal  model  used  in  [10] .  The 

phase  angles  of  a  were  assumed  to  be  independent  and  uniformly 

distributed.  To  evaluate  the  Cramer-Rao  bound  on  jointly 

estimating  t_  and  a,  one  is  required  to  evaluate  the  (3N  x  3N) 

Fisher  information  matrix,  [10]  . 

In  this  report,  we  assume  that  the  phase  difference 

between  two  successive  targets,  A0 .  .  ,  satisfies  the  following 

1  /  1  »  1 

relation. 


A0.  ... 
l ,  l+l 


(Ti+i-  V 


(2.4) 


where  A  is  the  wavelength.  With  the  above  assumption,  one  is  only 
required  to  jointly  estimate  the  magnitude  of  the  scattering 
amplitudes  | a | ,  and  the  phase  angle  of  the  first  target,  0^.  The 
associated  Fisher  information  matrix  has  dimension  { (2N+1) x  (2N+1) ) . 
It  may  be  written  as 
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Notice  that  we  have  neglected  the  details  of  the  lower  half  of^ 
since  it  is  a  symmetric  matrix.  The  inverse  of  &  is  the  Cramer- 
Rao  lower  bound  on  the  covariance  of  jointly  estimating  [  a  j  ,  t_, 
and  6 ..  . 

We  note  that  the  expectation,  E [ • ]  used  in  Equation  (2.5) 
applies  to  all  random  variables  of  the  likelihood  function  A . 

It  was  shown  in  [17]  that  a  tighter  bound  can  be  obtained  if  the 
expectations  in  Eq.  (2.5)  are  conditioned  on  the  random  parameters; 
and  the  expectation  with  respect  to  the  random  parameters  is  taken 
after  taking  the  inverse  of  the  Fisher  information  matrix.  For 
example,  the  phase  difference  of  two  scatterers  was  assumed  uni¬ 
formly  random  in  [10].  Let  (^"denote  the  Fisher  information  matrix 
obtained  with  the  expectation  taken  over  all  random  variables  and 
^A0  the  Fisher's  information  matrix  with  conditional  expectation 
upon  AO,  then 

EA0  [#0]  l2-6> 

where  E^g [  ]  denotes  the  expectation  with  respect  to  AO.  The 
above  property  was  used  in  [10].  We  will  discuss  its  application 
to  our  case  in  the  numerical  result  section. 
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3.  THE  CRAMER -RAO  BOUND  FOR  A  TWO-TARGET  MODEL 

In  this  section,  we  discuss  the  problem  of  the  two-target 
case.  The  expression  can  be  easily  extended  to  the  multiple-target 
case.  Rewriting  Eq.  (2.2)  explicitly  for  the  two-target  case  one 
obt„ ■ ns 


J  =  N  Jin  A 
o 


iO,  ,  ,  iA0.  ,  -i0i  .  * ,  „  *  -iA0. 

=  e  1  (a1Y1+cc2Y2e  )  +  e  1^aiYi  +a2Y2  e  * 

2  2 

-(a^  +a2  2a^a2p  cos  A0) 


(3.1) 


where  =  |  a^  |  , 

A©  "  02-0l  =  4(T2-T1) 

p(T2,T1)  =  Js  (t-X.^  s  (t-T2)dt 

Notice  that  for  most  radar  signals,  p(t2,t^)  =  p(t2~t^)  A  p(x). 

We  use  p  and  p  to  denote  the  first  and  second  derivatives  of  p 
with  respect  to  x,  respectively.  The  root-mean-square  (rms) 
bandwidth  (the  normalized  second  central  moment  of  the  signal 
spectrum)  of  s(t)  is  denoted  by  8  and  is  equal  to  the  square  root 
of  the  negative  of  p  evaluated  at  x=0.* 

With  the  above  definitions,  we  now  proceed  to  define 
the  Fisher  information  matrix. 


*For  a  linear  FM  signal,  8  is  approximately  equal  to  4ttB/c/12 
where  B  is  the  LFM  bandwidth  in  Hz,  and  c  is  the  speed  of  light. 


Z\  M 
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[~3  2  g,  nA  1  =  2_ 
[3t30iJ  No 


0  A  77  4  77  •  ~ 

—  £ ot 2  ( +  °‘1“2[  (r~)p  cos  A®“P  sin  aoi 


a22(^)  +  a^ct2^X~)P  cos  A0+p  sin  A0) 


(3.6) 


-E 


[3  2HnA~] 

bo,2  J 


2  22 

=  rr— (a,  +a_  +2a.a_  p  cos  A0 ) 
N  1  2  12 

o 


(3.7) 


Notice  that  the  above  expressions  do  not  depend  upon 
0 .  ,  the  initial  phase.  The  relative  phase  AO  is  a  function  of 

target  separation.  The  expectation  is  taken  with  respect  to  the 

4 

receiver  noise  only.  The  inverse  of  the  above  matrix  is  the 
lower  bound  on  the  covariance  of  jointly  estimating  a,  t_,  and  0^ 
For  the  purpose  of  comparison,  the  Fisher  information 
matrix  for  the  case  of  jointly  estimating  |  a  j  ,  i_,  and  0^  is  shown 
in  Appendix  C. 

It  is  always  interesting  to  know  the  case  when  the 
targets  are  disjoint  in  time,  i.e.,  p (t ) =p (t ) =p  (t ) =0 .  If  we 
were  to  assume  AO  random  the  information  matrix  would  be  diagonal 
(see  Appendix  C)  but  this  is  not  true  for  a  coherent  phase 
relation.  We  discuss  this  case  below. 
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3.2.  Results  for  Two  Resolved  Targets 


Substituting  p=p=p=0  into  equations  (3.2)— (3.7) ,  one 
finds  that  the  estimation  of  a^  and  are  uncorrelated  with  the 
estimation  of  t_  and  0^.  The  lower  bound  on  the  variance  of 
estimating  can  be  expressed  as 


,  2>!e 

“i  -  2 


a . 

l 

SNR. 

1 


(3.s; 


where  SNR.  = 

l 


2a  .  ‘ 

l 


N 


This  result  is  the  same  as  the  random  phase 

v-/ 

case.  The  covariance  of  estimating  t_  and  0^  is  lower  bounded  by 
the  inverse  of  the  following  matrix. 


32i>nA  i 

9  2£nA 

2 12  : 

1 

3t301 

—  p 

1 

1 

d2inh  ! 

3F^Tt 

3  ^InA 
3012 

Z  2q  2  ,  2.4TT.2 

al  6  +a2 

2  , 4  tt 
~a2  (  — 

)2 

2  , 4  7T ,  2 

-a2  (— ) 

a2262+a22 

<?>2 

2  ,4tt,  2 
“a2  ()T) 

2,4- 

a2  A 

) 

2  ,4tt 
2  (A  } 

n  2  ,  4tt  . 
a2  (A-) 

2  2 
°1  +a2 


(3.9) 


After  some  tedious  manipulation  one  finds  that  the 
bounds  on  the  variance  for  t ^ ,  t 2,  and  At  A  t2_ti  are 
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2.  (ot12+a22)a2232+a12a22(^)2 

°Ti  ~  2,  2.  2  20  2  fQ  2  ,  ,4ir .  2 . 

1  (a^  +a2  )a1  a2  $  [6  +  (— )  ] 


(3.8) 


.  2  2.  2„  2  ,  2  2  ,  4tt  .  2 

2  ,  *  1  +  2  ^  1  ^  +al  a2  *A  * 
°x0  -  ,  2  2,  2  20  2  r  D  2  ,  ,  4rr ,  2  , 

2  ( a  ^  2  a2  ^  ^  "^  ( i  )  ] 


lA 


0.9: 


oA  2  > 
At  — 


62(a12+a22)2 


(ai2+a22)ai2a22e2  [B2+  (^)  2) 


(3.10) 


Letting  R  A  — -  =  scatterer  amplitude  ratio  and 

1  2a.2 

SNR^  A  — —  =  peak  signal-to-noise  ratio  of  i-th  scatterer 
o 

one  can  rewrite  the  above  expressions  to  a  more  familiar  form. 


SNR^2  (1+R2) 


I1  +  -42. 


(3.8a) 


SNR282 (1+R2) 


R2  + 


-4>2j 


(3.9a) 


o.  >  — 
At  -  „2 


i+4>2,  t3™1  +  S™J 


(3.10a) 


Comparing  the  above  results  with  the  case  of  random  phase  (see 
Appendix  C) ,  one  notices  that  the  difference  is  in  the  appearance 
of  the  wavelength  dependent  term  (4tt/AB).  If  this  term  does  not 
appear,  Eqs.  ( 3 . 8a ) - ( 3 . 10a )  will  be  identical  to  the  case  when 
the  relative  phase  is  uniformly  random. 
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3.3.  Numerical  Results 


In  this  section,  we  present  some  numerical  results  and 
compare  them  with  those  for  random  phase.  We  assume  that  the 
signal  autocorrelation  function  is  represented  by  a  Gaussian 
pulse  shape 


P  (t)  =  e 


2  2 
-t  B  /  2 


Its  first  and  second  derivatives  are 


PIT)  =  -T82e-'2fi2/2 


2,2 


(x)  =  -B2  (l-x232)e*T  3  /2 


(3.11) 


(3.12) 

(3.13) 


In  order  to  facilitate  the  presentation  of  our  numeri¬ 
cal  results,  we  use  the  following  change  of  variables. 

(1)  The  Cramer-Rao  bound  will  be  evaluated  as  a  function 
of  target  separation  normalized  with  respect  to  the  inverse  of 
signal  root-mean-sqvare  bandwidth,  i.e., 


(2) 


i  A 


The 


xB 


(3.14' 

normalized  target  separation 

relative  target  phase  can  be  rewritten  in  terms  of 


i  as 


AO 


4  XT  T 

~r~ 


(3.15) 
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(3)  Examining  the  Fisher  information  matrix,  it  is 


evident  that  one  can  normalize  it  with  respect  to  0.  This  results 
in  changing  all  4tt/X  to  4tt/XB. 

(4)  The  final  results  will  be  shown  as  degradation  factors, 
i,  ,  the  Cramer-Rao  bound  with  closely  spaced  targets 
normalized  by  the  Cramer-Rao  bound  with  completely  resolved 
targets.  This  eliminates  the  explicit  dependence  on  the  signal 
bandwidth  and  wavelength,  rather,  the  results  can  be  evaluated 
using  the  product  of  wavelength  and  bandwidth  XB,  as  a  parameter. 

Generally,  the  signal  bandwidth  in  frequency  units  (e.g., 
the  Linear  Frequency  Modulated  (LFM)  bandwidth)  is  at  most  equal 
to  10%  of  the  center  frequency.  Let  this  bandwidth  be  denoted  by 


B  =  kf 


(3.16) 


where  k  is  a  positive  constant  less  than  or  equal  to  .1,  f  is  the 
center  frequency,  and  c  is  the  speed  of  light. 

The  above  bandwidth  is  related  to  the  root-mean-square 
bandwidth  6 (m  ^)  by  the  following  equation 


6  = 


c^  1 2 


c  /12 


(3. J  rt) 


l 


Equivalently,  one  has 

4  tt  __  2/1"  (3.19) 

A  5  k 

For  example,  if  one  considers  a  10%  bandwidth,  one  uses 
k  =  . 1.  This  results  in 

—  =  20/3  (3.20) 

a  r> 

A  5%  bandwidth  case  will  have 

yj  =  40/3  (3.21) 

The  degradation  factors  for  estimating  the  first  target 
amplitude,  location,  and  relative  target  separation  for  the  equal 
target  amplitude  case  with  5%  bandwidth  are  shown  in  Figs.  3.1- 
3.3,  respectively.  The  corresponding  results  for  the  10%  band¬ 
width  case  are  shown  respectively  in  Figs.  3. 4-3. 6.  Notice  that 
these  curves  show  a  general  decrease  with  increasing  target 
separation  with  a  superimposed  ripple.  Using  Eqs .  (  3 . 19 ) - ( 3 . 2 1 ) 

and  these  figures,  it  is  clear  that  the  ripple  period  for  5% 
bandwidth  is  exactly  naif  of  that  for  of  10%  bandwidth.  Also 
traced  are  the  results  for  the  random  phase  angle.  It  is  evident 
that  the  coherent  phase  information  is  able  to  substantially 
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reduce  the  estimation  error. 


Degradation  Factor  (dB) 


(Relative  Separation )  x  (RMS  Bandwidth ) 


Fiq.  3.2.  Degradation  factor  for  location  estimates,  5% 
bandwidth . 


Degradation  factor 


CQ 

X> 


Fig-  3.3. 


•5  .75 

•Relative  Separation) 


*•  1. 25 


1 5  1.  75 


ao 


bandwidth.  De9radati0n  factor  fQr 


xiRMS  Bandwidth) 

Oration  estimates. 


5% 


Degradation  Factor  (dB) 


(Relative  Separation)  x  (RMS  Bandwidth) 

Fig.  3.4.  Degradation  factor  for  amplitude  estimates,  10% 
bandwidth . 
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Degradation  Factor  (dB) 


.25  .5  .75  1.  1.25  1.5  1.75  2  0 

(Relative  Separation)  x  (RMS  Bandwidth) 


Fig.  3.5.  Degradation  factor  for  location  estimates,  10% 
bandwidth. 


•®  .75  1.  1.25  1.5  1.75  Z  0 


(Relative  Separation)  x  (RMS  Bandwidth) 

Fig.  3.6.  Degradation  factor  for  separation  estimates,  10% 
bandwidth . 


2: 


To  get  a  feeling  for  the  average  estimation  performance. 


we  may  treat  the  relative  phase  as  a  random  variable  uniformly 
distributed  between  0  and  2tt  for  the  bound  computation.  In  this 
case,  we  used  the  modified  bound  developed  in  [17]  to  take  the 
expectation  of  the  inverse  of  the  Fisher  information  matrix 
with  respect  to  A0.  This  step  also  eliminates  the  dependence 
on  the  percent  bandwidth.  The  results  for  amplitude,  location, 
and  separation  estimates  are  shown  in  Figs.  3. 7-3. 9,  respectively. 

There  may  be  situations  where  even  though  the  phase 
difference  is  coherent,  the  estimator  does  not  make  use  of 
this  knowledge.  The  Cramer-Rao  bound  for  this  situation  is 
obtained  by  treating  phase  angles  as  unknown  parameters  (Ref.  [10] 
or  Appendix  C)  and  not  averaging  over  them.  That  is,  the  bound 
of  Ref.  [10]  is  evaluated  for  phase  difference  as  a  function  of 
target  separation.  This  step  makes  the  bound  depend  upon  the 
wavelength  and  bandwidth  product.  Results  for  5%  and  10%  band- 
widths  are  shown  in  Figs.  3.10-3.15.  Notice  that  these  results 
are  extremely  oscillatory.  The  5%  bandwidth  case  has  a  period 
exactly  equal  to  half  of  10%  bandwidth.  These  results  should  be 
compared  with  the  curves  labelled  "with  random  phase"  of  Figs. 

3. 1-3. 6. 

Comparing  the  above  results,  one  concludes  that  using 
the  information  contained  in  the  phase  difference  can  substantially 
reduce  the  error  of  target  amplitude  and  location  estimation. 

The  maximum  likelihood  estimator  which  asymptotically  approaches 
these  performance  bounds  is  discussed  in  Appendix  B. 
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Fig.  3.7.  Degradation  factor  for  amplitude  estimates,  averaged 
over  relative  phase. 
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Fig.  3.8.  Degradation  factor  for  location  estimates,  averaged 
over  relative  phase. 
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Fig.  3.10.  Degradation  factor  for  amplitude  estimates  by 
neglecting  phase  information,  5%  bandwidth. 
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Fig.  3.11.  Degradation  factor  for  location  estimates  by 
neglecting  phase  information,  5%  bandwidth. 
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Fig.  3.12.  Degradation  factor  for  separation  estimates  by 
neglecting  phase  information,  5%  bandwidth. 
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Fig.  3.13.  Degradation  factor  for  amplitude  estimates  by 
neglecting  phase  information,  10%  bandwidth. 
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Fig.  3.14.  Degradation  factor  for  location  estimates  by 
neglecting  phase  information,  10%  bandwidth. 
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Fig.  3.15.  Degradation  factor  for  separation  estimates  by 
neglecting  phase  information,  10%  bandwidth. 
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4 .  SUMMARY 


In  this  report,  we  have  discussed  the  accuracy  of  parameter 
estimation  of  radar  signals  closely  spaced  in  range  using  the 
Cramer-Rao  bound  analysis.  Specifically,  we  extended  the  result 
of  [10]  to  include  the  case  when  the  phase  angle  from  successive 
signals  are  coherent.  Use  of  this  additional  information  is  shown 
to  substantially  reduce  the  estimation  error. 

A  comparison  of  the  performance  of  these  estimators  measured 
by  the  target  separation  required  to  achieve  a  3  dB  degradation 
is  summarized  in  Table  4.1.  Notice  that  because  the  degradation 
curves  are  oscillatory,  nominal  or  a  range  of  values  is  given 
for  some  cases. 

Although  there  are  no  simulation  results  presented  in  this 
report,  it  is  known  that  the  maximum  likelihood  estimator  can 
achieve  these  performance  bounds  when  the  signal-to-noise  ratio 
is  high,  [10],  [15],  [16],  For  this  reason,  the  maximum 
likelihood  estimator  equations  for  the  radar  signal  resolution 
problem  are  presented  in  Appendix  B. 
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TABLE  4.1 


SUMMARY  OF  PERFORMANCE  COMPARISON 


^SParameter 

Normalized*  Target  Separation  for  3  dB  Degradation 

Phase 

Relation  N. 

Amplitude 

Location 

Separation 

Coherent 

0.9 

—  0.35 

~  0.8 

Random 

1.45 

1.45 

1.45 

Coherent  but 
Ignored 

1.15-1.65 

1.2-1. 7 

0. 8-2.0 

*Separation  is  in  units  of  1/RMS  Bandwidth  or  equivalently 
0.55  x  nominal  range  resolution  (C/2B) .  Rayleigh  criterion 
corresponds  to  normalized  separation  of  1.8. 
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APPENDIX  A: 
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r  n 

ui  =  m|s(M.)  Y] 

j=l 


aj*s(t-v 


M 

dt+  z/  (t-fi)nm* (t) dt  =  M  g i 

m=l 


(A. 3) 


where  q 


i  £ 


L  j-i 


a . *s  (t-t  . ) 
3  3 


M 

dt+  *Z 

m=i 


Notice  that  the  q.  above  is  similar  to  the  y.  defined  in  (2.2). 

l  l 

Their  first  terms  are  the  same  and  their  second  terms  have  the 
same  variance  Nq/2 . 

Using  the  above  notation  in  (A. 2)  yields 


=  N  £nA 
o 


M  £aT3+a*\*-aTpa*  ] 


(A. 4) 


From  the  above  derivation  it  is  evident  that  the  Cramer-Rao 
bound  for  measurements  containing  m  pulses  is  equal  to  the  Cramer-Rao 
bound  for  a  single  pulse  divided  by  M. 

We  note  that  the  above  derivation  uses  the  assumption  that 
all  pulses  can  be  aligned  properly  for  the  required  summation. 

It  also  assumes  that  the  relative  position  and  orientation  of 
all  scatterers  are  preserved  for  M  pulses. 
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APPENDIX  B: 


The  Maximum  Likelihood  Estimator 
The  maximum  likelihood  estimates  are  those  parameter  values 

maximizing  the  log  likelihood  function.  The  log  likelihood 
function  for  our  case  is 

AnA  =  jjp(aTI  +  a*TY*  -  aTPa*)  (B.l) 

o 

In  the  following,  we  will  discuss  the  cases  of  targets  with 
random  phase  and  coherent  phase  individually.  In  the  last  sub¬ 
section,  we  will  extend  our  results  to  the  case  of  multiple 
measurements . 


B.l.  Targets  with  Random  Phase 

Ir  this  case,  one  first  maximizes  (B.l)  with  respect 
to  a  this  yields  the  estimate  of  a 

a  =  P-1^*  (B. 2) 


Substituting  a  into  (B.l)  yields 

~  1  T  -1 

w  =  rip  i*  (B. 

o 

The  t_  estimate  is  therefore  the  vector  ?  which  maximizes  (B.2). 
For  the  two  target  case,  the  above  expression  becomes 


inh  = 


[Y1Y1*+Y2Y2*‘P  ^2_tl)  (Y1Y2*+Y2Y1*)  1 


No(1-p  (W> 


(B .  4) 


The  estimates  (f^,^)  are  obtained  by  searching  through  the  set 
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of  twc -target  (after  matched  filtering)  time  samples  which 
maximizes  the  above  expression.  The  estimator  structure  is 
illustrated  in  Fig.  B.l. 

The  above  algorithm  assumes  that  the  received  signal  is 
analog.  If  only  sampled  data  is  available,  one  may  first  use 
the  available  samples  to  find  an  approximate  solution  and  refine 
this  solution  by  interpolation.  One  may  also  devise  algorithm 
working  directly  with  the  sampled  data.  One  such  algorithm  is 
defined  in  Ref.  [14],  We  will  not  discuss  this  further  here. 


B. 2 .  Targets  with  Coherent  Phase 

In  this  case,  one  does  not  have  to  estimate  all  complex 
target  amplitudes,  but,  rather  the  first  phase  angle  and  the  magnitude 
of  all  target  amplitudes. 

Let 


a  = 


“l 


N 


$  = 


iO  •y  -l 

e  <■  >  i 


0 

>ie  3 , 1 


0 

0 

0 


6 . eiGw'1 


where  0.  1  =0.-0,  =  ^-(x  . -t  . )  .  Then  a  =  $  a 

j/i  3  1  A  3  1  —  — 
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Fig.  B.l. 


Estimator  structure. 


The  .og  likelihood  ratio  becomes 

J  =  ei01aT<J>y  +  e~i0 1aT<t*y *  -  aT$P4>*a 
One  may  maximize  (B.5)  with  respect  to  a,  this  yields 
fi  =  [<5>P<i>*  +  4>*P$]~1[ei01<i>y  +  e”1014>*y*] 


(B.5) 


(B .  6 ) 


(B .  7 ) 


(B .  8 ) 


The  maximum  likelihood  estimates  of  are  those  values 

which  maximize  (B.8). 

For  the  two-target  case,  one  obtains  the  following 
expression  after  some  tedious  manipulations. 

„  [d  2+d  2  -  2d,d-p )cos  AG] 

J  =  — - - - - - - - -  (B .  9 ) 

N  (1-p^ (t--f . )cos  A© ) 

O  2  1 

where  A0  = 

A 

-  1  y x i  coste^+d x) 

a 

d2  -  |  y2I  cos  (G^tAOn-^) 

6^  =  phase  angle  of  y^ 


The  estimator  structure  for  the  above  problem  is  similar 
to  that  of  the  previous  subsection.  The  difference  is  that  one 
has  to  search  for  three  parameters  0^  t ^ ,  and  t2  instead  of  two. 


B. 3 .  With  Multiple  Pulses 

Similar  to  the  case  for  the  Cramer-Rao  bound,  the  log  like¬ 


lihood  ratio  becomes 


-  ■  ‘it  4  *  -At 


ifii 


*|  -  M*aTPa* 


(B.10) 


Maximizing  J  with  repect  to  a  yields 


=  p  1.4 


‘(s  §  -) 


(B.ll) 


Substituting  (B.ll)  into  (B.10)  yields 

;  ■  -t  -Ms  -) 

/\ 

The  optimal  estimate  ?  is  obtained  by  maximizing  J  with  respect 
to  t. 

Notice  that  the  above  estimator  first  requires  a 
summation  of  all  pulses.  This  assumes  that  all  pulses  can  be 
aligned  properly.  This  assumption  may  not  be  true  in  general 
due  to  the  performance  limit  of  the  realtime  tracking 
algorithm.  One  alternative  is  to  process  each  pulse  individually 
then  average  the  results  (e.g.,  for  the  target  separation  estimate). 
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APPENDIX  Cj 


The  Fisher  Information  Matrix  for  a  Two-Target  Model  With  Randan  Phase  Angles 
For  the  purpose  of  comparison,  we  state  the  Fisher  infor¬ 
mation  matrix  for  a  two- target  model  with  random  phase  angles 
in  this  appendix.  This  case  was  first  discussed  in  Ref.  [10]. 
Using  (2.3)  and  (3.1)  one  has 


j  =  N  £nA 
o 


_  „  v  +  a  y  +  a  v 

—  c*  ^  I  ^  2  ^  2  1 ' 1  /  ^ 

(C.l) 

-  a,2  -  a02  -  20(^02  cos  (0,-0  •,  )  P  ( T  2- T 1 ) 


The  Fisher  information  matrix  has  the  following  terms. 


'  2  1 

1  pcos  A0 

E 

2 

=  2 

6  a/ 

pcos  A0  1 

t 

r 

-J— 1  = 

PL2] 


~E[$ 


E[aa3x] 


-a^a2pcos  A0 


2  2 

-a-j^pcos  A0  a2  B 


a^^cos  A0  p (x ) 


=  2 


[ct^a 2cos  A0  p  (x ) 


a2cos  A0  p (x) 


cos  A0  p 
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Notice  that  if  two  targets  are  completely  resolved  so  that 
p(t)  =  p(t)  =  0,  then  the  above  matrix  becomes  a  diagonal  matrix 
and  one  obtains  the  following  familiar  equations 
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